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SUMMARY 


A comprehensive compilation of available turbulent-flow separation data for over- 
expanded supersonic nozzles is presented with a discussion of correlation techniques and 
prediction methods. Data are grouped by nozzle types: conical, contoured, and two- 

dimensional wedge. Correlation of conical-nozzle separation is found to be independent 
of nozzle divergence half-angle above about 9°, whereas the contoured-nozzle data follow 
a different correlation curve. Zero-pressure-gradient prediction techniques are shown 
to predict adequately the higher divergence-angle conical separation data, and an empiri- 
cal equation is given for the contoured-nozzle data correlation. Flow conditions for which 
the correlations are invalid are discussed and bounded. A nozzle boundary-layer transi- 
tion criterion is presented which can be used to show that much of the noncorrelating data 
in the literature are concerned with nonturbulent separation and which explains the pre- 
viously reported '’external flow effects” on nozzle separation. 

INTRODUCTION 

The flow field of a supersonic nozzle exhausting to an ambient pressure greater 
than that for which the nozzle is designed will contain oblique shock waves, and the flow 
may separate from the walls as a supersonic jet which fills only a portion of the cross- 
sectional area available for expansion. This realistic behavior is in contrast to the theo- 
retical inviscid description (ref. 1) of supersonic expansion after the throat, followed by 
a normal shock and subsonic compression. This is the classical flow pattern that would 
exist if there were no boundary layer. Figure 1 illustrates the difference between the 
two cases. 

For the case with separation, the nozzle flow expands from the supply pressure p^ 
to a minimum wall static pressure, usually designated in nozzle-separation studies as the 
separation pressure p g . Separation occurs upstream of where a normal shock would 
occur and, when the flow separates, the pressure rises to ambient pressure much more 
rapidly than the pressure rise associated with subsonic compression. This results in 
higher wall pressures for the separated case and improved performance (higher thrust) 
over the normal shock case. For vehicles such as proposed hypersonic research aircraft 
which have afterbodies that act as the nozzle, separation location will affect lift and stabil- 
ity characteristics as well as thrust. Obviously, knowledge of the point of separation is 
essential for performance prediction. Studies have been made (refs. 2 to 4) of techniques 
to induce and control separation for both improved stability and performance. 



Nozzles for launch vehicles are usually designed for some average expansion ratio 
which would give the best overall performance through the full trajectory. Thus, the 
nozzle will be operating at an off-design overexpanded condition at lower altitudes and 
will be susceptible to separation. 

The present paper presents a comprehensive review of available supersonic-nozzle 
separation data. A review of correlation and prediction techniques is made, and predic- 
tions are compared with the separation characteristics for various types of supersonic 
nozzles. It is shown that much of the scatter in separation pressure ratios at low Mach 
numbers is due to the existence of laminar or transitional flow which also explains the 
difference in data taken with and without an external coflowing stream. An appendix is 
included which gives an example of the application of the correlation techniques to a typi- 
cal nozzle. 


SYMBOLS 

A cross-sectional area of nozzle 

A* cross-sectional area of nozzle at throat 

f function 

M Mach number 

Mj Mach number ahead of separation shock 

M 2 Mach number behind separation shock 

p pressure 

R^- r transition Reynolds number 

R x Reynolds number based on length along nozzle from throat 

R* Reynolds number based on throat (sonic) conditions and throat diameter 

u s velocity at boundary- layer edge 

Ug characteristic velocity in boundary layer (ref. 31) 
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X 


spatial coordinate along nozzle axis with origin at throat 


a half-angle of conical nozzle or local wall angle for other nozzles (see fig. 1) 

y specific heat ratio 

Subscripts: 

a external ambient conditions at nozzle exit 

e exit 

ext conditions in external coflowing stream at nozzle exit 

s undisturbed value just ahead of separation (see fig. 1) 

t total conditions in stagnation chamber 

w wall conditions 

REVIEW OF SUPERSONIC-NOZZLE SEPARATION DATA 

History 

Separation in nozzles was first noticed in the early 1900’s as steam turbines came 
into use; however, no criteria were established for the prediction of the onset of separa- 
tion. Detailed investigations did not appear until the advent of rocket propulsion in the 
late 1940's when nozzle flow separation became a matter of practical importance. Some 
of the first quantitative correlations related separation in conical and two-dimensional 
nozzles to settling-chamber pressure pj. (ref. 5). Most of these studies showed that 
the separation pressure ratio (ratio of the nozzle wall pressure just upstream of separa- 
tion to the external ambient pressure p s /p a ) was independent of design expansion ratio, 
nozzle divergence angle, ratio of specific heats, and gas temperature; as a consequence, 
a value of p s /p a = 0.4 was used as a separation criteria and is still quoted today (e.g., 
see ref. 6) although more recent studies have shown it to be inadequate. The results of 
Scheller and Bierlein (ref. 7) conflicted with the other early investigations in that they 
found the separation pressure to be dependent on nozzle divergence angle. They also 
suggested that Reynolds number might be an important correlation parameter , but the 
lack of available data at that time prevented any correlation. Meleney and Kuhns (ref. 8) 
reported that as separation location approached the throat, the separation pressure ratio 
was more influenced by the ambient pressure. Considerable experimental work has been 
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devoted to nozzle separation. The reported investigations are too numerous to mention 
individually, but references 9 to 49 in addition to those previously mentioned are exam- 
ples of the work done to determine the effects of various parameters on separation in 
supersonic nozzles. 


Correlation Parameters 

A number of correlations have been suggested by various investigators. Early 
workers in the field used 



and 



( 2 ) 


Green (ref. 13) suggested a modification of equation (2): 

Pa £s_ f /£a\ 

Pt Pt \Pt/ 


(3) 


which through the use of the redundant parameter p a /p^ resulted in apparent suppression 
of scatter. However, this reduction in scatter was actually accomplished through a con- 
version to a higher reference value. 

More recent investigations have suggested 


and 



(4) 



(5) 


Although the first of these equations (eq. (4)) has the merit of being written in more con- 
venient terms for the designer, with p^. being a function of engine design and p a a 
function of flight environment, the second equation (eq. (5)) is more closely related to the 
local conditions at boundary-layer separation and is used in the present paper. Lawrence 
(ref. 42) achieved better correlation with equation (5) than with equation (4) and suggested 
it to be the preferable correlation. The two equations are actually related through the 
isentropic relation 
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(6) 



with y having only a small effect at lower Mach numbers where data are available for a 
range of y. It will later be shown that the effect of y is not apparent from the availa- 
ble data. 


Prediction Methods 

A number of predictions and correlations of the separation pressure ratio for tur- 
bulent boundary layers are shown in figure 2. Also shown is the pressure ratio across a 
normal shock wave. In general, shock-separated boundary layers are classified in two 
groups: free and restricted separations. Chapman, Kuehn, and Larson (ref. 50) refer 
to the first of these as a free-interaction separation and define it as a separation which is 
not directly influenced by downstream geometry. Separation in overexpanded nozzles of 
wide divergence is an example of free-interaction separation since there is no reattach- 
ment to the nozzle and no required separated-flow direction. 

Reshotko and Tucker (ref. 51) developed a semianalytical prediction method for the 
onset of turbulent separation. In this method, empirical incompressible separation cri- 
teria are used in conjunction with a compressibility transformation to predict separation 
for the compressible case. A relationship between the local Mach number and velocity 
profile (form factor) provides a means of determining a Mach number ratio across the 
discontinuity for shock-induced separation. The Mach number ratio was found to be 0.762 
for zero-pressure-gradient flows and somewhat lower for favorable pressure-gradient 
flows. Calculations of P s /P a as a function of M s are shown in figure 2 for this method 
using the Mach number ratio of 0.762 with y = 1.2 and 1.4. This technique and its appli- 
cation to separation in nozzles are discussed more fully in Lawrence's thesis (ref. 42). 

Arens and Spiegler (ref. 34) utilized the assumption, first suggested by Gadd 
(ref. 52), that the pressure rise associated with separation must be sufficient to stag- 
nate a characteristic velocity in the boundary layer u|. This theory is also shown in 
figure 2 for a value of u|/u s = 0.6 (as suggested in ref. 34) where u s is the local 
velocity at the edge of the boundary layer at separation. Reference 34 points out that 
this theory (as well as all others shown herein) assumes the separation peak pressure to 
be equal to the ambient exhaust pressure and does not account for any compression asso- 
ciated with the mixing region downstream of separation. 

Also shown in figure 2 are the empirical curves representing the zero-pressure- 
gradient flat-plate data of Chapman, Kuehn, and Larson (ref. 50) and Sterrett and Emery 
(ref. 53) for incipient separation. Reference 50 observed a slight Reynolds number effect 
on the separation pressure ratio, and the curve shown represents data for a Reynolds 
number of 10 6 . Reference 53 did not show an effect of Reynolds number and Holden 


5 



(ref. 54) showed the effect of Reynolds number to be on the order of R^VlO a j. higher 
Reynolds numbers (above 10 6 ). The pressure rise associated with a normal shock as 
well as the theory of Reshotko and Tucker (ref. 51) is shown (see fig. 2) for comparison 
on all subsequent correlation plots of p s /p a as a function of separation Mach number. 

CORRELATION OF SEPARATION DATA FOR SUPERSONIC NOZZLES 
Factors Causing Data Scatter 

Transition effects .- Some typical examples of nozzle- separation data are corre- 
lated in figure 3 in terms of P s /P a as a function of M s . The data of Meleney and 
Kuhns (ref. 8) show a wide variation in separation pressure ratio at low Mach numbers 
(see fig. 3(a)), typical of the scatter found in the literature. It is assumed by the present 
authors that the data for the low ambient pressure (0.33 atm or 33.43 kPa in fig. 3(a)) 
below M s = 2.4 are for laminar or transitional separation. The data of Herbert and 
Martlew (ref. 38) are reported to be laminar for the smooth-wall case and to be tripped 
by roughness ahead of the throat for the turbulent case (noted as rough-wall data in 
fig. 3(b)). 

In order to reduce the scatter in data correlations, especially at the low Mach num- 
bers, an attempt was made to define a boundary between laminar and turbulent separation 
by utilizing the data of Meleney and Kuhns (ref. 8). This was done by cross-plotting the 
data to allow the separation pressure ratio to be determined as a function of Reynolds 
number for constant values of separation Mach number M s . Two Reynolds numbers 
were calculated: the first was a length Reynolds number R x , based on free- stream 
conditions at the point of separation and axial distance from the throat; the second was 
a throat Reynolds number R*, based on sonic flow conditions and the throat diameter. 
Some authors have suggested R* as a characteristic of nozzles (ref. 37). 

The separation pressure ratios as a function of each of these Reynolds numbers are 
shown in figure 4. The data for each separation Mach number are characterized by a 
region of only small changes in separation pressure ratio until some lower value of 
Reynolds number is reached and the separation pressure ratio suddenly increases. 

This sudden break in the curve is believed to be the demarcation between laminar and 
transitional separation. The Reynolds number at this point of sudden increase is shown 
in figure 5 against separation Mach number and is defined as a transition Reynolds num- 
ber. Data falling below these curves are considered to be cases of laminar separation. 
Although this criterion is based on only one set of data, it does eliminate much of the 
scatter in the correlation of data as will be subsequently seen. This criterion is applied 
to nozzles of widely differing shape since there are not enough data to establish the effect 
of geometry on the transition Reynolds numbers. Reference 8 does not distinguish 
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between the data as laminar or turbulent but simply refers to the low Mach number var- 
iation in the separation pressure ratio as being "influenced to a greater extent by the 
discharge pressure p a ." By the reasoning given herein, it is felt that the supply pres- 
sure p^ is the primary influence since it determines the Reynolds number. 

Nozzle-exit effects .- The data of figure 3(b) show another feature of nozzle sepa- 
ration which has increased data scatter in proposed correlations. Separation near the 
nozzle exit causes the separation pressure ratio to rise above the normal downward trend 
with increasing M s . This region of increasing P s /P a is small for the turbulent data 
and somewhat larger for the laminar data. In an effort to correlate the location in the 
nozzle at which the ratio P s /P a increases above the expected trend, the nozzle pres- 
sure ratio P^/P s and separation Mach number at which "exit deviation" occurs from the 
correlation of P s /P a against M s are shown as a function of nozzle-design pressure 
ratio Pt/Pe aiK * design exit Mach number in figure 6. There is scatter in the data but 
the curve given by 



often suggested in the literature (ref. 40) seems to be a reasonable guide above a Mach 
number of 2. A modification of equation (7) given as 

= 0.8 - l\ + 1 (8) 

A* \A* j 

seems to fit the data better at low exit Mach numbers. Data falling below this curve 
should follow the correlation of P s /P a plotted against M s . 

Also indicated in figure 3(b) are example values for which an oblique shock exists 
at the nozzle exit without separation occurring (M s = M e ). The limiting overall nozzle 
pressure ratio Pt/P a for this condition is shown in figure 7 as a function of design pres- 
sure ratio. These data can be well represented by the curve (ref. 19) given by 

= 1 + 0.39 (9) 

Pa p e 

Pressure ratios p t /p a above this curve will not separate the flow. 

Correlations of Experimental Data 

The available values of separation pressure ratios for nozzles are shown in fig- 
ures 8 to 11 as a function of separation Mach number. Experimental data are included 
for conical and contoured, two-dimensional, and axisymmetric nozzles exhausting into 
still air and with a coflowing external stream. Data which showed the exit deviation 
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noted earlier in the discussion of figures 3(b) and 6 are omitted from figures 8 to 11. 

Data for which Reynolds numbers could be calculated and which fell below the transition 
Reynolds number curves of figure 5 are shown as solid symbols. For much of the availa- 
ble data (see table I) there is not sufficient information to calculate the Reynolds numbers. 

Conical nozzles .- Separation pressure ratios for conical nozzles with low divergence 
angles shown in figure 8(a) are seen to be in fair agreement with the theory of Reshotko 
and Tucker (ref. 51) only around a Mach number of 2 and to fall considerably below the 
theory at higher separation Mach numbers. The data of figure 8(b) for slightly higher 
values of nozzle divergence more closely approach the theory throughout the separation 
Mach number range above M s = 1.7. The turbulent data for all higher divergence angles 
(figs. 8(c) and (d)) closely agree with the theory above M s = 1.7. Below M s = 1.7 the 
data for all divergence angles fall below the theory, thus showing a tendency to level off 
at a value of P s /P a between 0.5 and 0.6. There may be a stronger Reynolds number 
effect in this region where, in general, the Reynolds numbers are very low. Figure 4 
shows some variation of separation pressure ratio with Reynolds number even for the 
turbulent data. This deviation from the theory may also be a result of the very strong 
pressure gradient close to the throat of the nozzle. As was indicated earlier, a lower 
value of the Mach number ratio used in the Reshotko and Tucker theory (ref. 51) was 
suggested for favorable pressure-gradient flows. Mager (ref. 55) also indicates that the 
ratio should decrease slightly at lower Mach numbers. 

Note that a few data points show a pressure rise at separation greater than that for 
a normal shock loss. (In figs. 8(a) and (b), p s /p a is smaller than normal shock val- 
ues.) This occurs because the entire pressure rise is not obtained at the point of sepa- 
ration. If a normal shock exists at the separation point for the data below the normal- 
shock pressure-rise curve, the flow then diffuses subsonically to the ambient exhaust 
pressure. (See fig. 1(b).) 

The data of figures 8(b) to (d) have specific heat ratios of 1.2 to 1.4 (see table I), 
but no effect of specific heat ratio on separation pressure ratio is discernible in the data. 
The theories of Arens and Spiegler (ref. 31) and of Reshotko and Tucker (ref. 51) show an 
increasing effect of specific heat ratio on the separation pressure ratio as Mach number 
increases (fig. 2). The available data having a specific heat ratio of other than 1.4 are 
limited to data below M s = 3.2, and a specific-heat-ratio effect may be masked by the 
data scatter. 

Two-dimensional wedge nozzles .- Separation pressure ratios for two-dimensional 
wedge-flow nozzles (fig. 9) correlate similarly as those for conical nozzles. The cause 
of the rise in the separation pressure ratio above M s = 3.1 for McKenney's data (ref. 11) 
is not known but it is not a low Reynolds number effect. Much of the scatter in the data 
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below M s = 2 may be due to low turbulent Reynolds numbers or laminar flow; however, 
the Reynolds numbers are not available for the data of reference 31. 

Axisymmetrical contoured nozzles .- The data for contoured nozzles (fig. 10), the 
case of most interest for modern nozzle design, fall below the prediction of Reshotko and 
Tucker (ref. 51) at separation Mach numbers greater than about 3. A similar trend was 
found with low-divergence-angle conical nozzles (fig. 8(a)) which is not surprising since 
both types of nozzles have low local wall angles above M s = 3 (near the exit for con- 
toured nozzles). The close proximity of the wall may eliminate these cases from the 
free- interaction separation category discussed earlier. Lawrence (ref. 42) and Guman 
(ref. 56) suggest that the close proximity of the wall to the shear layer (see fig. 1(a)) 
causes a change in the entrainment process and that the pressure rise is the pressure 
necessary to drive the reverse-flow entrainment air in the narrow region between the 
wall and the mixing region. Lawrence also showed that using the pressure plateau just 
downstream of separation in place of p a resulted in agreement with theory. Guman 
(ref. 56) suggested the use of a percentage of p a (65 percent) as a correlating parameter 
but this would probably vary with each case. A knowledge of the entrainment process is 
necessary for a prediction technique to be applied. 

An empirical equation is shown in figure 10 which is a second-order curve fit to the 
data (other than ref. 33) between M s = 2.4 and M s =4.5. This equation is 

^ = 1.082 - 0.363M S + 0.0386M g 2 (10) 

The data of Roschke and Massier (ref. 33) deviate rapidly upward above M s = 3.4 
and are never in close agreement with the rest of the data. Although the Reynolds num- 
bers are not available for these data, it is noted in reference 33 that the rapidly rising 
data were obtained at low back pressures, which indicates lower supply pressures p^. 
for the same overall pressure ratios. Thus, the flow may have been laminar at these 
lower supply pressures. (See fig. 3(a).) 

Nozzles with external coflowing stream .- Several studies have been made of nozzles 
exhausting into coflowing external flows. The first of these (ref. 15) indicates that for 
supersonic external streams the flow separates at much lower back pressures (higher 
p g /P a ) than with quiescent external air. These data are shown in figure 11 along with 
the data of references 22, 27, and 37, which would agree with the conclusion of refer- 
ence 15 until the Reynolds numbers are examined. Elimination of the laminar points 
and the data of reference 15, for which Reynolds numbers are not available, gives results 
which follow a trend similar to that of the higher divergence-angle conical nozzles. 
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Comparison of data for various nozzles .- A summary of the data of figures 8 to 11 
is shown in figure 12. Here it can be clearly seen that the correlation of Reshotko and 
Tucker (ref. 51) is valid for the high-divergence-angle conical nozzles and is invalid for 
the low-divergence-angle conical nozzles and contoured nozzles at higher separation 
Mach numbers. Also shown is the empirical equation (eq. (10)) for the contoured noz- 
zles. The appendix gives an example of the application of the correlation techniques to 
a typical nozzle. 


CONCLUSIONS 

An analysis of available experimental data on turbulent-flow separation in over- 
expanded supersonic nozzles resulted in the following conclusions: 

1. The data can be correlated for each type of nozzle if data near the nozzle exit 
are excluded and if laminar data are distinguished and eliminated. 

2. Separation in nozzles can be predicted by using zero-pressure-gradient free- 
interaction theory over most of the nozzle length for wall divergence angles greater than 
about 10°. 

3. Prediction of separation for low-divergence-angle conical nozzles and for the 
low wall-angle region of contoured nozzles requires knowledge of the pressure rise in the 
entrainment process and is no longer a free-interaction separation. 

4. No obvious effect of specific heat ratio was found when the data were correlated 
in the form of separation pressure ratio as a function of separation Mach number. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
May 31, 1978 
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APPENDIX 


EXAMPLE OF USE OF CORRELATION 

A sample of the procedure for determining the separation point in a particular noz- 
zle at certain flow conditions is presented in this appendix. Figure A1 is a sketch of a 
typical conical nozzle having a 15° half -angle and an exit-to-throat area ratio of 8.125. 
The linear dimensions of the nozzle are not needed in applying the separation criterion 
if the flow is known to be turbulent. If not known to be turbulent, the length or throat 
Reynolds numbers may be calculated for comparison with the transition criterion of 
figure 5. 

Figure A2 gives the calculated pressure distribution for the nozzle if assuming one- 
dimensional flow and a specific heat ratio equal to 1.4. By using the theory of Reshotko 
and Tucker (ref. 51) which was a good fit to the 15° conical-nozzle data of figure 8(c) 
(other empirical fairings of data or theories could be used), the values of P s /P a a re 
determined for various Mach numbers along the nozzle at various x-positions. The sep- 
aration Mach number is simply the wall Mach number immediately ahead of separation 
and is a function of the nozzle design. The value of P s /P a a t each x-position divided 
by the corresponding value of p w /p t (fig. A2) gives the desired ratio of nozzle total 
pressure (a design parameter) to ambient pressure (a function of flight environment). 

The resulting curve of separation location as a function of Pj./p a is shown in figure A3. 
The position downstream of which the correlation no longer holds (from fig. 6) is indi- 
cated near the exit as point A, and the value of p t /p a (from fig. 7) above which separa- 
tion will not occur is indicated at the exit as point B. A straight-line fairing between the 
two values is shown but the actual curve may differ somewhat. 
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(a) Separation case. Viscid flow. 

Pa-"^ 


(b) Normal-shock case. Inviscid flow. 


Wall 

pressure, 

p s 

x 

(c) Wall pressure distributions. 
Figure 1.- Flow in overexpanded nozzles. 
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Figure 6.- Pressure ratio at which turbulent separation pressure ratio P s /P a 
deviates from correlation. Mach number is calculated for y = 1.4. For 
symbol key see table I. 




Figure 7.- Minimum nozzle pressure ratio Pj/P a for which the oblique shock 
remains at exit and flow does not separate. Turbulent flow. For symbol 
key see table I. 




Figure 8.- Separation pressure ratio as a function of separation Mach number. For symbol key see table I 

Solid symbols denote laminar data as defined in figure 5. 






Figure 8.- Continued. 
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Figure 9.- Separation pressure ratio as a function of separation Mach number for two-dimensional 

1 o n 

wedge-flow nozzles, a = 2± to 30 . For symbol key see table I. 




Figure 10.- Separation pressure ratio as a function of separation Mach number for axisymmetrical 
contoured nozzles. For symbol key see table I. Solid symbols denote laminar data as defined 
in figure 5. 



Figure 11.- Separation pressure ratio as a function of separation Mach number for tests 
with a coflowing external stream. Solid symbols denote laminar data as defined in 
figure 5. 
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Fijfui 6 12,- Summary plot of separation pressure ratios as a function 
of separation Mach number. 
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